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In the mean field approximation of (2+l)-flavor Nambu-Jona-Lasinio model, we strictly derive 
several sets of coupled equations for the chiral susceptibility, the quark number susceptibility, etc. 
at finite temperature and quark chemical potential. The critical exponents of these susceptibilities 
in the vicinity of the QCD critical end point (CEP) are presented in SU(2) and SU(3) cases, respec¬ 
tively. It is found that these various susceptibilities share almost the same critical behavior near the 
CEP. The comparisons between the critical exponents for the order parameters and the theoretical 
predictions are also included. 
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I. INTRODUCTION 

It is believed that the strongly-interacting matter un¬ 
dergoes a phase transition to the quark-gluon plasma 
(QGP) at high temperature T and/or at high quark 
chemical potential fi. Partially restored chiral symme¬ 
try and deconfined quarks and gluons are expected in 
this new state of matter [1-6]. The phase transition 
should have implications for the study of the on-going 
heavy-ion collision experiments at the BNL Relativistic 
Heavy-Ion Collision (RHIC) and the Large Hadron Col¬ 
lider (LHC) [7-10] and the future FAIR(GSI) project in 
Darmstadt and NIC A project in Dubna [11-13]. 

The phase transition from the hadronic matter to the 
QGP at finite T with vanishing fi has been investigated 
extensively on the lattice and within the effective the¬ 
oretical models. It is believed that with two massless 
quarks the chiral phase transition is of second order and 
the critical behavior falls into the universality class of 
the 0(4) Heisenberg magnet spin model in three dimen¬ 
sions [14, 15]. However, in nature, the u- and d-quark 
have small but finite masses and the second order phase 
transition turns into a smooth crossover. 

The reliable lattice simulations of the QCD chiral tran¬ 
sition with finite /i are not available due to the severe 
fermion sign problem. Many effective theories [1, 16- 
19] of QCD predict the existence of the tricritical point 
(TCP) and CEP in the QCD T — {i phase diagram in 
massless and massive quark cases, respectively. The TCP 
is the point where a line of critical points (0(4) line, a 
second-order transition line) at lower /i’s and higher T’s 
turns into a first-order phase transition line at higher jiA 
and lower T’s. With the current quark masses increas¬ 
ing from zero to the physical masses, the TCP emerges 
to the CEP at which the first-order phase transition line 
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ends from higher /i’s towards lower /i’s in the phase dia¬ 
gram. The location of the CEP has been reported to be 
accessible at RHIC [20]. 

Both at the TCP and CEP, second-order phase tran¬ 
sitions are considered to occur, which are characterized 
by the long-wavelength fluctuations of the order param¬ 
eter. Therefore, the values of critical exponents will be 
an interesting and important question, which is related 
to the universality class of this phase transition. Based 
on the analysis of the Landau-Ginzburg thermodynamic 
potential’s expansion in the order parameter, the crit¬ 
ical exponents in mean field approximation have been 
predicted [21]. Furthermore it’s expected that the crit¬ 
ical exponents of the susceptibility depends on whether 
the path to approach the GEP is asymptotically parallel 
to the first order transition line or not [22]. It is nec¬ 
essary to verify these universality arguments from the 
calculations of specific effective theory of QCD [23]. The 
Nambu-Jona-Lasinio (NJL) model is believed to be a 
good choice [19, 24-27]. 

The rest of this paper is organized as follows. In Sec. 
H the formulation of the (2+l)-flavor NJL model and the 
QCD thermodynamics are briefly reviewed. In Sec. HI 
we present the derivation of several sets of coupled equa¬ 
tions for chiral susceptibility, quark number susceptibil¬ 
ity, etc. at finite temperature and quark chemical poten¬ 
tial. In Sec. IV, the numerical results of several suscep¬ 
tibilities and their corresponding critical exponents are 
presented. It will be found that they are consistent with 
the predictions by the formalism of Landau-Ginzburg 
thermodynamic potential’s expansion in the order pa¬ 
rameter [21]. The discussions on the critical exponents 
for the order parameters at the CEP and vanishing chem¬ 
ical potential in the chiral limit are also included. Finally, 
in Sec.V we will summarize our results and give the con¬ 
clusions. 
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II. THE NAMBU-JONA-LASINIO MODEL {{q.Q)) quark number density {{q^q)) are given by 


The Lagrangian of the SU(3) NJL model [16, 28-31] is 
given by 

^ - m)ip + Gi[{ipXi'ip)‘^ + (ipiy5Xi'ip)'^] 

- G2[{i’Xi'j^'tpf + {'tpXa^j5'tpf] (1) 

- K{det[ipi{l + 75)V’i] + det[?^i(l - 75)V’j]}> 

where the column vector ^ = {u, d, s) is the quark field 
with three flavors, Nf = 3, and three colors, Nc = 3. 
m = diag(m^i, rrid^ rris) is the current quark mass matrix. 

are the Gell-Mann matrices, a=l,2,3,...,8, 

The second and third terms represent four-point scalar- 
pseudoscalar (SP) and vector-axialvector (VA) interac¬ 
tions, respectively, which are invariant under the chiral 
U(3)l G U(3)r symmetry. The last term represents six- 
point interaction which breaks the axial U(1)a symme¬ 
try, but it is invariant under SU(3)l G SU(3)r. These 
four-point and six-point interactions are constructed in 
accordance with all relevant symmetries of QCD. 

The color singlet part of the six-point interaction can 
be written 


^6 =^Dijk[^{ipXiip){ipXjip){ipXkip) 
- {■>Ph5Xitp){ipi^5Xjip){tpXkip)]. 


The totally symmetric coefficients Dijk are identical 
to the SU(3) structure constants dijk for i^j^k G 
{1,2,3,...,8}; Dojk = -l/Vd for j = k = 1,2,3,...,8 


and Dooo = y §• 

Developing the thermodynamics of ^njl iu the mean 
field approximation, we obtain the thermodynamical po¬ 
tential density uj = Q/V. The conditions duj/d{Sm) = 
duj/d{5ii) = 0, which minimize the thermodynamical po¬ 
tential, lead to the gap equations for each of the quark 
flavors i\ 


5mi = Mi-mi = -dGi{{qiqi)) + 2K{{qjqj)){{qkqk)), (3) 
with i^j^k cyclic, and 

5iii = iiri-IJ-i=dG2{{qlqi)). (4) 

Where Mi and firi are the constituent quark mass and 
renormalized chemical potential of flavor i, respectively. 
Here the thermal expectation value of an operator 0 is 
denoted as 


((©)) 


Tr 0 


(5) 


Where is the Hamiltonian operator, while is the 
quark number density operator. The chiral condensate 


{{qq}) = - n{p, Iir) - m{p, Iir)}dp, (6) 

N 

{{(lU)) = -^l P^inip^Pr)-m{p,iir)]dp, (7) 

Jo 

with 


n{p,p} 


1 

1 + exp [I3{E - p)] ’ 


(8) 


m{p, p) 


1 

1 + exp [/?(£' + p)] ’ 


(9) 


E = \J + p^ and /3 = 1/T. A is a momentum cut- 
off, which is introduced to regulate the ultraviolet di¬ 
vergence. For convenience we will use the non-covariant 
three-momentum cutoff scheme [32]. 

The thermodynamical potential density u = VtjV is 
given by 


=— log[Tr 

E j\\E,+T{\og[l 




d] 


i=u,d,s 


-h log[l-h e'^^ ^''>])}dp + 2Gi{{qiqi)f 

- ‘^G2{{qlqi)f - 4:K{{ququ)){{qdqd)){{qsqs))- 


( 10 ) 


In this paper, the widely accepted parameter set will 
be employed [33]: mu,d = 5.5 MeV, rus = 140.7 MeV, 
A = 602.3 MeV, Gi’a^ = 1.835, G 2 = 0, KA^ = 
12.36, which yields = 92.4 MeV, = 135.0 MeV, 
Mk = 497.7 MeV, = 514.8 MeV, M^/ = 957.8 MeV 
and quark condensates {{duqu))^^^ = {{ddqd))^^^ = 

—241.9 MeV, {{qsqs))^^^ = —257.7 MeV. The constituent 
quark masses are also obtained: Mu = M^ = 367.7 MeV 
and Ms = 549.5 MeV. Here we don’t include the vector- 
axialvector (VA) interaction. However, for the general¬ 
ization of the equations to calculate the susceptibilities 
shown below, we keep the G 2 coupling in our following 
derivation. 

Now Eqs. (3),(4),(6) and (7) form a set of self- 
consistent equations. By solving these equations itera¬ 
tively, one can obtain the effective quark masses, chiral 
condensates and quark number densities at finite tem¬ 
perature and quark chemical potential (only the simple 
physical case jdu = jdd = jdg is considered in this pa¬ 
per). In the co-existence region of Wigner solution (chiral 
symmetry partially restored) and Nambu solution (chi¬ 
ral symmetry broken), the solution which minimizes the 
thermodynamical potential density uj is stable [25]. The 
thermodynamical potential density uo = ft/V is given in 
Eq. (10). 

Having taken this into account, we show the final re¬ 
sults in Eig. 1 and Eig. 2. It’s worth noting that, at 
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FIG. 1. (color online): Up/Down quark mass 
[MeV] 



FIG. 2. (color online): Strange quark mass M. 
[MeV] 


T=0, Ms stays almost constant after /j,crit where first- 
order transition occurs until /i exceeds Ms and we start 
to have density of strange quarks in the system. However, 
at finite T, we start to have density of strange quarks at 
a relatively small /i in the system and this is a differ¬ 
ent feature compared with the T=0 case. According to 
our numerical calculation, the number density of s-quark 
is approximately one tenth of that of u(d)-quark near 
the CEP. Therefore we can expect that the presence of 
s-quark will influence the location of the CEP and the 
interplay of s-quark with u(d)-quarks has the potential 
to affect the critical behavior near the CEP. Hence we 
study critical behavior related with the strange sector in 
the vicinity of the CEP in the following. 


III. EXACT EQUATIONS FOR VARIOUS 
SUSCEPTIBILITIES 

Experimentally, the linear responses, such as the sus¬ 
ceptibility, conductivity, etc., of the physical system to 


some external field, are often measured to study the prop¬ 
erties of the related system. Therefore the studies of 
various susceptibilities, as linear responses of the quark 
matter to the external fields are very important on the 
theoretical side, which characterize the non-perturbative 
properties of QCD vacuum [34-38] and are also widely 
used to study the phase transitions of strongly interact¬ 
ing matter [39-42]. 

Now let us introduce the definitions of six kinds of sus¬ 
ceptibilities to be discussed in this work: the chiral sus¬ 
ceptibility Xs, the quark number susceptibility Xq^ fhe 
vector-scalar susceptibility Xvs^ fhe thermal susceptibil¬ 
ity Xt and two auxiliary susceptibilities Xm and Xn- For 
mathematical convenience we first introduce four of these 
susceptibilities in the free quark gas case (the interac¬ 
tion terms in the Lagrangian are set to be zero, i. e., 
^int = 0) [43], where M and fir are reduced to m and 
li^ which are independent quantities. Denoted with the 
superscript their definitions and expressions are as 
follows: 


^( 0 ) = _ dim))} 

Xs — 


dm 


Jo 




dp, 




(0) _ ^((gg))/ 

Jo 


Xvs — 


y(0) = 

Am 


dji 

djqU))} 

dm 


7r2 , 

^ ./o 


E 


-h{ii)dp, 


y(0) 

A.VS f 


( 11 ) 

( 12 ) 

(13) 

(14) 


where ^(p) + ft.(p) = 2n(p)(l - n(p)), ^(p) - ft.(p) = 
2m(p)(l — m.(p)), /(p) = 1 — ri(p) — m(p), the subscript 
“/” represents the free quark gas systems. It should be 
noted that x^m and have the same analytical expres¬ 
sion, which is reasonable from the viewpoint of statistical 
mechanics: 


Xm Xvs vdmdp ^ ’ 

where is the QCD partition function in the free quark 
gas case. 

In the interacting case. Mi and jiri {i = u^d^s) are 
no longer independent, they both are the functions of 
mi = mu = md, m2 = ms, Mi = Mn = Md, M2 = Ms 
and T. In the following we make the convention that 
iqiqi)) = imqu)) = iwa)) and iq 2 q 2 )) = iqsqs))- Here 
the SU(2) flavor symmetry has been required. Let us 
specify the notations of these susceptibilities, taking the 
chiral susceptibilities as examples, x^j = where 

i,j G {1,2}. These susceptibilities are coupled with each 
other as follows: 
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xV + 4Gix“ - {{q2q2))xV “ {{qiqi))xf) + xii°H4G2x“), (16) 

xf =xf °\4.Gixf - 4ii'((gigi))x“) + X?l°H4G2Xm), (17) 

Xm =xT\^ + 4Gix“ - ^K{{q2q2))xV “ 27sr((9i9i))xf) + xJ^°n4G2Xm), (18) 

Xm =xTK^Gixf - 4 iX((gigi))x“) + x^^°^( 4 G 2 Xm)- ( 19 ) 

xT =Xs^°ni + 4G'iXs^ “ 4K {{qiqi))xl‘^) + X«i°H4G2Xm)i (20) 

xf =x\^^\^Gixf - ‘^K{{qiqi))xf - 2K{{q2q2))xf) + xlT{^G2Xm)^ (21) 

X™ =X^°ni + 4Gixf - 4.K{{qiq{))xf) + xrn4G2X™), (22) 

Xm =Xm°H4Gixf - 2K{{qiqi))xf - 2K{{q2q2))xf) + xT\^G2Xm)- (23) 

xll =Xs^°^(4Gix“ - 2i^((92?2))x“ - {{qiqi))xll) + xil°^(l + 4G2x“), (24) 

xll =xT\^Gixll - iK{{qiqi}}xll) + xT\4:G2xf), (25) 

x“ =xTK^Gixll - ‘^K{{q 2 q 2 ))xll - {{qiqi))xll) + x\^^\'^ + 4G'2x“), (26) 

xf =Xm°^(4GiX?^ - 4iX((gigi))x“) + xT\'^G2xf)- (27) 

X«s =Xs^°H4GiX^s - 4iX((gigi))x),s) + X«1°^(1 +4G2 Xq^)) (28) 

xll =Xs^°n4Gixis - ‘^K{{qiqi))xll “ {{q2q2))xll) + xSH4:G2xl^), (29) 

Xq^ =Xm°K^Gixll — 4:K {{qiqi))xll) + Xq^°Hl +4G2Xq^)) (30) 

Xq^ =xlli°H4:Gixll - 2K{{qiqi))xll “ ‘^K{{q 2 q 2 ))xll) + Xq^°n4G2Xq^)- (31) 


XT^l^°\^GiXT-2K{{q2q2))XT-‘^K{{qiqi))x%)+xlT{^G2xll)+Mii3x\^°'>-Hrii3Xm^\ (32) 

Xt =Xs^°n4GiXT - ^K{{qiqi))xT) + X?l°n4G2X^) + M2/3Xq^°^ - Mr2/3x^°\ (33) 

x)i =x),^°n4GiXT-27f((g2g2))XT-2i^((9i9i))Xr) + Xq^°H4G2Xn)-Mri/3Xq^°^+xJ^°\ (34) 

Xn =Xm°H4G'iXr - 4iX((gigi))xT) + Xq^°H4G'2X^) - Mr2/3Xq^°^ + xf • (35) 


where Xr — Xn — The superscripts 

2 ( 0 ) represent the free susceptibilities of the u(d)-quark 
sector and s-quark sector, respectively. 

x(0) = ^ f p^E/3^h{iJ,)dp. (36) 

Jo 

The inputs of variables of the free susceptibilities denoted 
with superscript and functions mentioned in Eqs. 
(16)-(35) are constituent quark mass M and renormal¬ 
ized chemical potential (ir instead of current quark mass 
m and bare chemical potential ji. Now five sets of lin¬ 
ear equations for these susceptibilities are presented here. 
One can solve these sets of equations with the simple 
methods of linear algebra (set fii = /i 2 again), and then 
these susceptibilities will be at hand. The same as in the 
SU(2) case [44], we use the peaks of various susceptibili¬ 
ties xV^iXW Xt criteria to locate the critical 


crossover region. Fig. 3 reveals an obvious critical band 
in the crossover region. Hence the same conclusion as in 
the SU(2) case can be drawn in the SU(3) case: it is more 
suitable to define a critical band rather than an exclusive 
line in the crossover region. 


IV. CRITICAL EXPONENTS FOR VARIOUS 
SUSCEPTIBILITIES 


As is well known that a second-order phase transi¬ 
tion occurs at the CEP, where the correlation length 
tends to infinity and these susceptibilities diverge [19, 
21, 23, 24]. Here in the SU(3) case the CEP is located 
at (T^^P,/i^^P) = (67.73 MeV, 318.42 MeV) [24] with 
iququ)X^ = -191.9 MeV, = -251.8 MeV, 

liqUj) = 3.171 X 10® MeV® and (iqU,)) = 2.424 x 
10® MeV®, while in the SU(2) case, = 





5 


+ First-order transition 
☆ CEP 
o Peak of Xg^ 

X Peak of x"'^ 

Ays 

□ Peak of xl 


II [MeV] 


FIG. 3. (color online): Phase diagram with a critical band in 
the crossover region. 


(31.95 MeV, 346.64 MeV) with (iququY^^ = -206.5 MeV 
and (lYlqu)) = 2.431 x 10^ MeV^ using the parameter set 
and equations in the Ref. [44]. The critical behavior of 
the quantities of interest in the vicinity of the CEP can 
be well described by a set of indices, the so-called criti¬ 
cal exponents. The motivation for this study arises from 
the fundamental and universal phase transition consid¬ 
erations. Some analysis based on nonperturbative renor¬ 
malization group method makes very big yet not enough 
progress in the recent years [23, 45-52]. 

We adopt the quantity A((®gi)) = {{mi))w - {{mi))N 
or A((g|gj)) = {{q\qi))w - {{q\qi))N, along the first-order 
transition line towards the CEP, as the order parame¬ 
ter [18] by analogy with the liquid-gas phase transition, 
which falls into the Z 2 universality class, namely, the 3- 
dimensional Ising universality class [53, 54]. The sub¬ 
script N and W denote the Nambu phase (chiral sym¬ 
metry broken) and the Wigner phase (chiral symmetry 
partially restored), respectively. 

The determination of these critical exponents, which 
govern the strength of these susceptibilities’ divergence, 
depends on whether the path to approach the CEP is 
asymptotically parallel to the first order transition line 
or not [22]. Indeed, the temperature t and magnetic field 
H of the 3D Ising model are mapped as certain linear 
combinations of T and fi on the T — {i plane near the 
CEP. To be specific, the t direction is tangential to the 
first-order transition line while the H direction is not. 
If the path is parallel to the /i-axis in the T — /i plane, 
from lower (higher) /i’s towards the critical at fixed 

temperature T = (marked by ^ (^) in the follow¬ 

ing tables), the critical exponents (e') are defined as 
follows: 


Inxi = -ei In I yu - \ +Ci, (37) 

InXi = -e^\n \ y- y^^^ \ -hq. (38) 

where q is a constant hereafter. 

We can also choose the path parallel to the T-axis to 
approach the CEP and define the corresponding critical 


exponents following the same convention as above (the 
direction from lower T’s towards is marked by t 

and, conversely, ^ ). 

According to the definition of the baryon number sus¬ 
ceptibility, in the SU(3) case. 



E( 


d{{qUi)) 

dyi 


)t = 


X 


22 

<7 


3 


(39) 


Since x]^ is at least 2 orders of magnitude larger than 
near the CEP in our calculation, the critical behavior of 
Xb is mainly governed by While in the SU(2) case, 

Xb = x“- 

The critical exponents (3i for the order parameters are 
defined as follows: 


In A((5,«)) = |i„ KT - + (m - I +C,, 

(40) 

In A((,;i)) = ^In I (T - f | +c,, 

(41) 


where the path to approach the CEP is along the first 
order transition line. 

Eor comparison with the 0(4) spin model universality 
class predictions, we also consider the critical behavior 
of the {{qiqi)) {i = u^d) in the chiral limit at vanishing 
chemical potential and from lower T’s towards the critical 
temperature Tc (located at 193.22 MeV and 191.80 MeV 
in SU(2) and SU(3) case, respectively) : 


In I {{qiqi)) |= /?»In | T - Tc | +Ci. (42) 


The critical exponents for the order parameters both 
in SU(2) and SU(3) cases are presented in Table 1. The /3 
values of A {{qiqi)) and A{{qlq-)) towards the CEP along 
the first-order phase transition line are consistent with 
the theoretical estimated /3 value ^ in mean field theory 
instead of the p value 0.30 — 0.34 for the 3D Ising model 
[54]. Logarithmic fitting procedures are shown in Eig. 4 
taking the SU(3) case as an example. The [3 value of 
iiduqu)) at the critical temperature Tc and the vanishing 
chemical potential in the chiral limit is 0.50 (0.49) ±0.01 
in the SU(2) (SU(3)) case, which is also consistent with 
the predicted value ^ in mean field theory instead of the 
0(4) spin model. These results are reasonable intuitively 
since the mean field approximation has been applied to 
the NJL model in this work. 

The critical exponents of the various susceptibilities in 
the SU(3) case for the u(d)-quark sector are presented 
in Table II, III and those for the s-quark sector in Table 
IV and partly in Eig. 5. The critical exponents in the 
SU(2) case are also presented in Table II, III and partly 
in Eig. 6. It can be seen from Table II-IV that these 
various susceptibilities share almost the same critical be¬ 
havior near the CEP, that is to say, the e values of these 
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TABLE L Critical exponents for order parameter 


Quantity 

Critical exponents 

SU(2) NJL 

SU(3) NJL 

Universality 

A{{qnqu)} 

pu 

0.49±0.02 

0.51±0.01 

1/2 


Ps 


0.50±0.02 

1/2 

Miqtqu)) 

Pi 

0.48±0.02 

0.51±0.02 

1/2 

Miqlqa)) 

Pi 


0.50±0.02 

1/2 

^ququ)) 

Pu 

0.50±0.01 

0.49±0.01 

1/2 


TABLE IL Critical exponents (parallel to fi axis) 


Quantity 

Critical exponents/Path 

SU(2) NJL 

SU(3) NJL 

Xq' 

^ 

0.67±0.01 

0.69±0.01 


4^7 ^ 

0.68±0.01 

0.69±0.01 


4^ / ^ 

0.69±0.01 

0.70±0.01 



0.66±0.01 

0.65±0.02 

xl\ 

4s / ^ 

0.68±0.01 

0.71±0.01 


4s7 ^ 

0.67±0.01 

0.69±0.01 

Xt 

Ct 1 ^ 

0.65±0.01 

0.67±0.02 


47^ 

0.69±0.01 

0.71±0.02 

Xn 

4 / ^ 

0.64±0.01 

0.65±0.02 


47^ 

0.69±0.02 

0.70±0.01 


TABLE III. Critical exponents (parallel to T axis) 


Quantity 

Critical exponents/Path 

SU(2) NJL 

SU(3) NJL 



0.68±0.01 

0.70±0.01 



0.67±0.01 

0.68±0.01 


e^/t 

0.69±0.01 

0.71±0.01 


erv; 

0.66±0.01 

0.64±0.01 

xll 

e“/t 

0.68±0.01 

0.70±0.01 


e“7; 

0.67±0.01 

0.68±0.01 

Xt 

/ t 

0.66±0.01 

0.68±0.01 


/ 4- 

0.68±0.01 

0.70±0.01 

Xn 

ej. /t 

0.66±0.01 

0.66±0.02 


n 

0.69±0.01 

0.69±0.01 


TABLE IV. Critical exponents (for s quark sector) 


Quantity 

Critical exponents 

SU(3) NJL(|| /X axis) 

SU(3) NJL(|| T axis) 

xf 

4^ 

0.68±0.04A 

0.68±0.04/t 


22/ 

0.65±0.02/^ 

0.65±0.02/; 

xf 

22 

0.67±0.03/-s- 

0.67±0.03/t 


22/ 

^s 

0.66±0.02/^ 

0.65±0.02/; 

xll 

22 

^vs 

0.67±0.03/-s- 

0.68±0.03/t 


22/ 

^vs 

0.65±0.02/^ 

0.65±0.02/; 

Xt 


0.69±0.02/-s. 

0.70±0.02/t 


€t 

0.70±0.01/^ 

o.70±o.oi/; 

Xn 

'-n 

0.68±0.02/-^ 

0.69±0.03/t 


'-n 

0.68±0.01/^ 

0.67±0.02/; 


quantities are very close to the mean-field exponent | 
no matter whether the path is parallel to the T-axis or 
/i-axis and what the direction is chosen to approach the 
CEP both in SU(2) and SU(3) cases. It is reasonable 
since the critical fluctuation mixes in the quark number 
density fluctuation and chiral condensate fluctuation as 
well as in the thermal fluctuation at the CEP, due to fi¬ 
nite current quark mass m and quark chemical potential 
ji [21, 55]. This mixture can be well understood from 


Eqs. (16)-(35) by the coupling relations of and 

Xt (corresponding to the quark number density, chiral 
condensate and the thermal fluctuation, respectively), 
especially between different quark flavors. In Ref. [56] 
the path dependencies of these critical exponents near 
the CEP are also discussed in (Polyakov-)Quark-Meson 
model, which both confirm our result, and the interest¬ 
ing scaling crossing phenomenon in a narrow regime near 
CEP is demonstrated there. It’s worth pointing out that 
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FIG. 4. (color online): Order parameters A{{qiqi)), A((g|g^)) 
{i = u, s) as a function of (T^^^ — T)^ + _ ^^2 along 

the first order transition line towards the CEP in the SU(3) 
case. 



FIG. 5. (color online): Susceptibilities for the s-quark sector 
as a function of — T at in the SU(3) case. 



FIG. 6. (color online): Susceptibilities for the u(d)-quark 
sector as a function of — /x at in the SU(2) case. 

the critical exponents of the baryon number susceptibil¬ 
ity obtained in Ref. [19] are reproduced here and they 


agree with the value of the mean field theory prediction 
I [21]. In addition, Fig. 5 shows a crossover of different 
universality classes from the changes of the exponents 
related with the s-quark sector, which is consistent with 
the spirit of Ref. [21]. Therefore, it can be indicated that 
s-quark may has a smaller critical region for the CEP 
than u/d-quark does, from the comparison with the ex¬ 
ponents related with the u(d)-quark sector presented in 
Fig. 6. 


V. SUMMARY AND CONCLUSIONS 

In this paper we have strictly derived several sets of 
coupled equations for the various susceptibilities at fi¬ 
nite temperature and quark chemical potential in the 
mean field approximation of (2+l)-flavor Nambu-Jona- 
Lasinio model for the first time. Then the rationality of 
using these susceptibilities as the criteria to determine 
the crossover region in the SU(3) case is discussed. It 
is thought to be more suitable to define a critical band 
rather than an exclusive line in the crossover region. 

We have located the CEP at = 

(67.73 MeV, 318.42 MeV) for Nf = 3 and 

= (31.95 MeV, 346.64 MeV) for Nf = 2. 
The various susceptibilities and their corresponding 
critical exponents have been calculated in the SU(2) 
and SU(3) cases, respectively. For the first time our 
calculations roughly verify the mean-field prediction 
that near the CEP these susceptibilities share the same 
critical behavior no matter whether the path is parallel 
to the T-axis or /x-axis and what the direction is chosen 
to approach the CEP, which is due to the mixture of the 
quark number density fluctuation and chiral condensate 
fluctuation as well as the thermal fluctuation at CEP in 
the presence of finite current quark mass m and quark 
chemical potential /x. The critical exponents of the 
order parameters towards the CEP along the first order 
transition line and towards the Tc at vanishing quark 
chemical potential in the chiral limit both agree with 
the mean-field predictions. More insightful work beyond 
mean field approximation calls for more applications 
of the renormalization group method to this issue 
within the NJL model or other nonperturbative effective 
theories. 
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